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1 Introduction 



The boundary value problem for linear operators in non-trivial regions leads 
to complications and quite often it is necessary to resort to the numerical 
analysis even for the simplest operators which possess a known kernel in the 
whole space, such as e.g. Laplace operator in the R n M. 

The possibility of a new approach to these problems appeared along with 
developing the theory of point interactions in quantum mechanics, firstly 
stimulated by the famous Kronig- Penney model ||, and systematically in- 
vestigated in || where self-adjoint extensions for Hamiltonian with point-like 
interactions have been constructed so that the explicit form for resolvent has 
been obtained for some physically significant systems. It is in reference || 
that it has been already pointed out that the limiting case of infinitely strong 
point interaction allows effectively to split the space in two separated regions 
that lead to two boundary problems on half line. 

This important trick has been successfully developed in ||, |5|, |6| and it al- 
lows one to write down the explicit expressions for the Green function of 
Schrodinger equation for the particle in one dimensional and radial boxes 
with Dirichlet and Neumann boundary conditions provided the appropriate 
problem in the whole space has been solved. The technique that has been 
successfully used for such a derivation is a direct summation of perturbation 
series (Dyson series J7J), effectively leading to geometrical progression due 
to the specific form of the perturbation. We would like to stress here, that 
numerous analytical results obtained up to now are related to quantum one- 
dimensional problems (as e.g., Krein's formula 0), effective one-dimensional 
problems after variables' separation ('radial problem'), and point-like inter- 
actions (see || U and recent book || for detailed references which include 
both solvable 5-perturbation cases and boundary value problems). 

The question naturally arising here is whether it is possible to generalize 
these constructions to higher dimensional boundary problems. As we will 
see, and this will be the main aim of this work, it can be done at least 
for Dirichlet boundary conditions for arbitrary linear evolution operators in 
topologically trivial region (homeomorphic to a ball in R d ), when assuming 
some natural conditions for the propagator of "free particle" to be fulfilled. 
It is worthwhile to mention here that we intend neither to construct the 
self-ajoint extension of appropriate singular perturbed operators H, nor to 
perform a detailed investigation of convergence properties of appropriate per- 
turbation series (as is well known the answer may be negative even for regular 
perturbations, see e.g. 0). To the contrary, we intend to propose explicit 
construction of Green function for non-separable case and demonstrate its 
validity by known examples. Generally speaking the problem we will solve 
can be formulated on an abstract level of the theory of linear operators by 
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introducing some generalization of projector operators, but we expect that 
it could only shadow simple foundation of the approach we use. Moreover, 
although, admittedly, the manipulations with perturbation series and the 
subsequent limit of infinitely large coupling constant are rather formal, we 
do not know another way to obtain the Green function representation for 
general boundary problem that we will construct in this paper. 



2 Series summation for singular perturbed 
system 

Let us have a linear evolution equation in R d of the form 

(|-£)*(x)=0, (1) 

with explicitly time independent operator £ acting on the function defined 
over R d and obeying Dirichlet boundary condition \l/(x)|r = 0, where T = 
{x : P(x) = 0} is the boundary (hyper-surface) of the region B under con- 
sideration, in which we seek the solutions of equation ([I]). 

We start from a consideration of the "free particle" , omitting the boundary 
condition, and assuming that propagator for that case, given by 

K°{x',x";t) =< xVV > 9(t), (2) 

is already known. Here 9 is the Heaviside unit-step function, incorporated 
into (0) to ensure causality property K(x',x";t) = 0, when t < 0. The 
propagator K possesses the composition property 

K(x', x"; t) = J dx 1 K(x!, x x ; t 1 )K(x 1 , x"; t - h). (3) 

R d 

which follows from the semi-group property with respect to the time evolu- 
tion, which in turn is automatically fulfilled for an explicitly time-independent 
operator £ and the decomposition of unity in an appropriate functional space, 
namely 

Id = £ \x) (x\ , (4) 

x 

where summation (integration) is performed over discrete (continuous) index 
enumerating states (see e.g, 0). These properties are natural for most of 
physically significant models so that our consideration is not very restrictive. 

Now we will emulate the boundary condition by introducing into ([!]) the 
additional singular potential term of the form V s = — 7<Jp(x). The general- 
ized 5-function used here is a distribution concentrated on the hyper-surface 
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P [[UJ. In the limiting case of 7 — > 00 the corresponding one dimensional 
problem turns out to be a Dirichlet boundary problem on the half line |§ 
(see the discussion in introduction). We will demonstrate explicitly that the 
same situation is met in higher dimensions. 

We will use the method expounded in || f|, |, ||, [Ll|], performing a perturba- 
tion expansion, starting from the formula for the propagator of the singular 
perturbed problem 

K(x',xV) =< x'V^V > 9(t). (5) 



The formal perturbation series over powers of V s can be constructed as in 
the quantum mechanics [12 and reads 



00 „ „ 

K s (x', x"; t) = K°(x', x"; t) + £ 7 n / dt x j dxiK°(x.', x l5 h - 0)5 P (xi) x 

n=l n r>H 



Rd 



n 

i=2 



J dtj J dXjK°(xj-i, Xj-; tj — tj_i)5 P {x.j 



Rd 



K°(x n ,x";t-t n ). (6) 



The convergence questions appearing at this moment should be treated for 
every problem considered separately, e.g. for Schrodinger equation the exis- 
tence of well defined Green function has been proven rigorously in some cases 
[[J. For arbitrary linear evolution equation we must stay on the formal level 
only to go further. 

After performing the Laplace transformation for the Green function, defined 
by 



G(x', x"; E) = J e~ Et K(x', x"; t)dt, 


the following series representation can be written 

00 . 

G 5 (x', x"; E) = G°(x', x"; E) + Y,l n d^G ^', xi; E)S P (x 1 

n=l 



(7) 



n 

3=2 



Rd 



G°(x n ,x";£). (8) 



The behaviour of G°(x',x"; E) at coincident space arguments in spaces with 
d > 2 may lead to divergence of integrals in (D), but we will not discuss this 
in details, since the appropriate procedures of regularizations are well known 
(for point-like perturbations in quantum mechanics see e.g., || |TT|). We 
only point out that for most interesting cases of two and three-dimensional 
quantum problems there are no singularities within our approach, opposite 



4 



to the models described in 01. Indeed, the short distance behaviour of the 



Green function in d-dimensional spaces is (|§, f.6.2.1.2) 

G(x',x",£) oc |x / -x"| 1 - d/2 yi- d /2(fc|x'-x"|), 



(9) 



where Y n (x) is Bessel function | 13fl , so that the relevant underlying singular- 
ities are integrable for d=2,3. For higher dimension (or) and other operator 
£, some sort of regularization should be used as e.g. in Of. 



Returning to our problem, now we can introduce new coordinates by the 
map F : x = {xi, ..x d } i— ► y = {y l7 ..,yd} with Jacobian J = so that 

the equation of hyper-surface P will be given by yd = r\ (see e.g.,|frOf). We 
designate all coordinates except the last one namely {yi : i = 1, ..d — 1}, by 
Q, so that y = {Q,y d }. Then, the integrations over ^-functions are simply 
projections on the submanifold, defined by yd = r\ and we get 

oo „ 

G 5 (x', x"; E) = G°(x', x"; £) + £ 7 ™ J y/g^G ^', (y d )', n u V ; E) x 



71=1 



n 

J'=2 



j y/gffiGPin^ri^w E) 



G (n n , V ,n",(y d y;E),(W) 



where the integration is performed over the hyper-surface P, g = det^"), 
g^ u = §^r|p is an induced metric tensor on P and we introduce coordinates 
Q, yd corresponding to the initial and final points x', x". Now we want to ex- 
pand the Green function G°(£lj-i, rj, Qj, r\\ E) in a series of functions defined 
on P. Let us choose an appropriate full (complete) orthonormal system of 
functions {f K (Q)} on the boundary T of the region B, where k is some multi- 
index enumerating the system /, with a standard L 2 (f2) scalar product 



(11) 



n 



For example, the case of axially symmetric closed surfaces has been recently 
treated by Prodan |TJ|], where the projection of the resolvent operator on 
such surfaces has been investigated. 



We represent G° in the form 

G°(n',tn",mE) = £ G K ',Atv;E)f K/ (n')f K/l (n'' 



(12) 



so that the coefficient G K > jK "(£, rj; E) is expressed as 

g K ,,A^mE) = I Jg(n') g (n») G°(cr, c, n", m E)f K ,(n')f K „{n") dn'dn" 



(13) 
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Then, substituting (|12D into ([T0|) we get 

oo - 

G 5 (x',x";E) = GV,x";£) + E7" / v^i^i E V, E) 

n=l p k',Ki 
Kj—\,Kj 



X 



/^(nOAiWII 

i=2 



L P 



G°(x', x"; £) + 7 E &Vu((lte)', 17; E)g Kn>K „( V , (y d )"; E)f K ,{ti)f K „{tt') 



x 



j fr\, ~\_ ^ h\, ft 



S K1 , Kn + lG KUKn {v, m E) + 7 2 E Gki, K2 (v, V, E)G K2 , Kn (r], v;E) + ... 



(14) 



where we used the orthonormality of the functions ([TT|). After summing up 
the geometrical progression we obtain 

G s (x!,x";E) = G°(x',x";E) + 

k' ,k" 

For the brevity of notation we used matrix form within the square brackets, 
and ( n /G) n is an ordinary matrix power. Taking the limit 7 — > 00 we get 
finally 

G 5 (x!,x"]E) = G°(x',x";E)- 
E [g((y d y^;E)g(r ] ,r ] ;E)- 1 g(r ] ,(y d )";E)} , „ f K '(tf)f K »(tf'), (16) 

K ,K 

which is the main result of our paper. 



3 Discussion 



As it is easily seen, the last formula solves an appropriate Dirichlet boundary 
problem. Indeed, the statement that the object constructed above satisfies 
the differential equation (|]) is evident from the construction and, if (yd)' = r\ 
or (yd)" = T], that is if the initial or the final points are on a boundary, the 
term in the square brackets simply gives the free Green function expansion 
coefficient and the whole sum becomes the free Green function cancelling 
the first term in (pi]), which means that the Dirichlet boundary condition is 
obeyed. 

It is worthwhile to point out here that the spectrum of the system under 
consideretion is given by such values of E that g(rj,r)]E) is non-invertible. 
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In two dimensional case the multi-index k becomes an ordinary one and we 
obtain the condition of vanishing determinant 



Det Q(rj, r\\ E) = 0. (17) 

From the last equation it is easily seen that our approach looks like some al- 
ternative and generalization of the boundary integral method JTR 16fl , where 
the spectrum of 2-D billiard can be obtained based on an integral of Green 
function's normal derivative over the boundary. 

It is also easy to demonstrate that formula ( ]T6|) leads to a known one for the 
case of the separability of variables. Let us do it explicitly for the 2-D case 
of quantum particle in a circular region (T = {x : |x| = R}). An appropriate 
formula for the Green function in polar coordinates (r, 0) reads (see also 
[ |T7|1 for alternative derivation) 

oo 

G(x,x';£)= J2 GfrS^eW-*), (18) 

m=— oo 

G (r r>- E) - C°(r r>- E) ^(r' , R- E)G^Ry- E) 

d(r, r,E)-G l (r,r,E) G o [R)R . E) ■ ( 19 ) 

A natural choice of the functions' family is of course f m {(p) = exp{im(j)}. 
Expanding the free particle Green function in the same manner as in (|18|) 
and calculating the coefficient of equation (0) one can see that 

g mm ,(r,r';E) = G m (r,r>; E)6 mm >, (20) 

so that the matrix inversion becomes trivial and after the substitution of (p0|) 
into (0), formula (|T9|) follows immediately. Similar arguments may be used 
for other separable quantum problems. 

Thus, we see, that indeed we have successfully constructed the explicit repre- 
sentation for the Green function of linear evolutional equation with Dirichlet 
boundary condition, based on Green function in the whole space, thereby 
generalizing results already known in separable cases in quantum mechanics. 

The formula fll6| ) can be rewritten in a more formal way, introducing the 
series expansion of (^(x^x"; E) in a manner like (12) for G°. Then, 



Gi>A(y d )'i (yd)"; e) = g K >A(y d Y, (y d )"; e) - 

g((y d )', rj; E)G( V , V ; E)" 1 ^, (y d )"; E)] , _„ (21) 



or in operator notation 



g s ((y d y,(y d y-,E) = g((y d y,( yd y-E)- 

g((y d )', 77; E)g(rj, m E)- 1 ^, (y d )"; E). (22) 



7 



The last expression is suitable for further formal manipulations in the case 
of double 5-perturbation V = j(8(y d — a) + 8(y d — b)), where the system is 
being "squeezed" into a narrow shell a < rj < b, simulating the quantization 
on a hyper-surface in the limit a — > b in a manner similar to [Q (see the 
equation 2.15 in [|J), but the detailed analysis will be published elsewhere. 

It should be mentioned also, that our approach can be modified also to be 
used for more general perturbation of the form V s = — 7/i(x)<5p(x) with 
arbitrary function h. Then similar arguments show that the only difference 
from the case considered above is that one should change yfg to y/g h(Q, rj). 
Then, we should use another function family for the expansion or, it may 
be more convenient to expand into the series not the Green function but the 
product 



h(n;OHW',v)G°(n',{i,n", m E) = £ g K/ ,Atv;E)fA^')fA^"), 

k" 

(23) 

and the final formula becomes 



(5V,x";£) = G°(x',x";£)- £ 



kW > yJh(W, (y d )')h(W, (y d )") 
g((y d )', rj; E)g( V , m EAGiri, (y d )"; E)] , „ . (24) 



Such a generalization may be useful for the systems with the boundary whose 
initial shape is not very convenient for the construction of the function fam- 
ily set f K and when it is easy to perform some transformations before using 
the proposed approach. In this case, after a transformation to the new co- 
ordinates (and, e.g., accompanied by the "local time rescaling" fL8|), the 
initially pure 5-function perturbation really transforms to a non-uniform one 
like discussed above. 
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